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ABSTRACT 


This  document  describes  and  documents  a  numerical  procedure  to  compile  data  sets  to 
be  applied  in  the  North  Atlantic  (NA)  basin-wide  numerical  simulations.  The  specific 
objectives  are; 

1)  to  acquire  climatological  data  sets  for  the  NA,  and  archive  them  on  a  data  library. 

2)  to  develop  an  interpolation  procedure  for  2-D,  and  3-D  data  fields  from  a  regular  grid 
into  the  generalized  coordinate  system  adopted  by  the  ocean  models. 

3)  to  prepare  plotting  facilities  for  the  display  and  presentation  of  the  original  and  inter¬ 
polated  data  sets. 

The  numerical  progreims  described  by  this  document  and  the  user’s  manual  associated  with 
them  were  developed  and  implemented  by  Ms.  Grace  Foo,  a  PhD  student  in  the  Scientific 
Computing  Program  at  USM,  and  part-time  research  assistant  with  COAM. 
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1.  INTRODUCTION 


Under  ONR  Grant  No.  N00014-92-J-4112,  The  Center  for  Ocean  Si  Atmospheric  Model¬ 
ing  (COAM)  at  The  Univer5ity  of  Southern  Mississippi  (USM)  has  undertaken  a  scientific 
progreun  to  develop  am  ocean  prediction  system  with  emphasis  in  the  North  Atlantic  that 
connects  a  basin-wide,  eddy-resolving  model  with  regional  models  incorporating  data  as¬ 
similation  procedures  to  provide  the  requested  forecasting  ability  in  selected  areas.  An 
essential  paxt  of  the  project  is  also  to  evaluate  new  physiceJ/numerical  formulations  that 
take  advantage  of  the  advancing  computer  technologies  and  are  suitable  for  basin-wide, 
eddy-resolving  ocean  circulation  simulations. 

Two  are  the  models  that  COAM  will  evaluate  and  apply  to  the  prediction  system;  The 
S-Coordinate  Rutgers  University  Model  (SCRUM),  and  the  Spectral  Finite  Element  (SFE) 
model,  both  designed  and  developed  by  Prof.  D.  Haidvogel  and  his  colleagues  at  Rutgers 
University. 

To  construct  a  robust  operational  prediction  system,  it  is  necessary  to  perform  extensive 
sensitivity  and  calibration  studies  to  evaluate  the  numerical  efficiency  and  physical  accu¬ 
racy  of  the  models.  Moreover,  both  SCRUM  and  SFE,  being  recently  released,  are  lacking 
a  background  of  testing  under  different  dynamical  regimes.  Therefore,  there  is  a  need  to 
design  a  series  of  preliminary  experiments  in  which  different  physical/numerical  issues  can 
be  isolated  easily  and  investigated  separately. 

This  document  is  focused  on  the  initialization  for  such  preliminary  simulations  and  de¬ 
scribes  the  interpolation  procedure  from  climatological  data  into  SCRUM  grid  configura¬ 
tions.  The  initialization  procedure  for  the  SFE  model  is  under  implementation.  Climato¬ 
logical  data  are  usually  on  a  coarse  but  regular  grid  that  removes  most  of  the  mesoscale 
signals  and  destroys  small-scale  features;  SCRUM  us»>s  an  orthogonal  curvilinear  coordi¬ 
nate  system  which  conforms  to  the  irregular  coastal  boimdaries  and  a  generalized  vertical 
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coordinate  trajisformation  that  follows  the  bottom  topography.  Numericed  experiments 
will  be  necessary  for  generating  a  numerical  grid  suitable  for  the  eissessment  of  the  basin¬ 
wide,  eddy-resolving  model.  Thus,  the  issue  is  to  implement  an  interpolation  package  that 
can  be  applied  easily,  is  numerically  efficient,  yet  preserves  the  dominant  features  of  the 
original  data  sets. 

This  document  describes  and  documents  a  numerical  procedure  to  compile  data  sets  to 
be  applied  in  the  North  Atlantic  (NA)  basin- wide  numericzJ  simulations.  The  specific 
objectives  are: 

1)  to  acquire  climatological  data  sets  for  the  NA,  and  archive  them  on  a  data  library 
(NADL). 

2)  to  develop  am  interpolation  procedure  for  2-D  and  3-D  data  fields  from  a  regular  grid 
into  the  generalized  coordinate  system  of  the  SCRUM  and  SFE  grids.  This  group  of 
programs  also  constitutes  the  initial  kernel  for  the  development  of  a  library  of  utility 
programs  (NAUL)to  be  used  in  connection  with  the  North  Atlantic  ocean  prediction 
system. 

3)  to  prepare  plotting  facilities  for  the  display  and  presentation  of  the  original  and  inter¬ 
polated  data  sets. 
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2.  NORTH  ATLANTIC  DATA  LIBRARY  (NADL) 


The  Levitus  annual  mean  of  temperature,  T,  and  salinity,  S,  and  the  Hellerman  wind 
stress  data  sets  constitute  the  first  kernel  for  the  North  Atlaintic  Data  Library  (NADL). 
The  library  will  contain  data  files  covering  the  region  extending  from  20  “S  (to  include  the 
Tropical  Convergence  Zone  (TCZ))  to  TS^N  (the  Fram  Straits)  from  iiO^W  to  20‘’E  (from 
American  to  European/ African  coasts)  (Fig.l). 

The  extracted  data  sets  are  usually  reformatted  according  to  the  stemdeirds  defined  for 
the  NADL.  The  NADL  is  intended  for  the  sole  internal  use  of  the  group  engaged  in  the 
development  of  the  ocean  prediction  syestem.  It  has  not  been  conceived  as  a  facility  to  be 
accessed  by  a  wide  number  of  users. 

3.  THE  INTERPOLATION  PROCEDURES 

The  Hellerman  and  Levitus  data  sets  eure  representative  of  2-D  vectorial  and  3-D  scalar  data 
fields,  respectively.  The  different  nature  of  the  original  data  sets  leads  to  different  issues 
zmd  approaches  in  the  interpolation  algorithms,  as  discussed  in  the  following  sections.  The 
results  are  presented  for  the  SCRUM  grid  illustrated  in  Fig.  1.  However,  the  interpolation 
procedure  is  independent  of  the  grid  configuration.  Information  relative  to  the  model 
domain  axe  extracted  from  the  files  generated  by  the  grid  generation  peickage  developed 
by  Hedstrom  and  Wilkins  (User’s  manual,  in  prepeu-ation)  and  supported  by  Prof.  D. 
Haidvogel  as  part  of  the  SCRUM  and  SFE  models.  Thus,  the  algorithms  described  in  this 
document  can  be  applied  also  to  the  SFE  model  with  miriimum  changes  in  the  input  file 
requirements. 


3 


3.1.  The  2-D  Vector  Field  Interpolation  Procedure 


Hellerman  winds  are  interpolated  on  two  steps  as  follows: 

1)  perform  a  2-D  horizontal  scalar  interpolation  of  the  (x.y)  (i.e,  longitude,  latitude)  wind 
components.  The  interpolation  is  done  using  the  bilinear  algorithm  developed  by  Dr. 
D.-S.  Ko  (Aug.  1989),  and  included  in  the  North  Atlantic  Utility  Library  (NAUL). 

2)  perform  a  rotation  of  the  axis  such  that  at  each  grid  point  the  wind  components  axe 
referred  to  the  orthogonal,  boundary  fitted  coordinate  system  of  SCRUM. 

The  annual  mean  Hellerman  wind  in  the  subdomain  stored  in  the  NADL  and  the  associated 
interpolated  field  are  depicted  in  Figs.  2.  and  3.,  respectively.  As  Figs.  2.  and  3.  indicate, 
the  program  changes  the  wind  stress  units.  The  Hellerman  wind  stresses,  (Tx,Ty),  are 
referred  on  CGS  units  (dyne/cm^);  on  the  other  hand,  SCRUM  uses  MKSA  imits  and 
takes  (tx,  ry)/po  .  Thus,  a  scale  factor  is  introduced  to  make  the  data  consistent  with  the 
model  requirements. 

Five  programs  are  used  to  process/plot  the  wind  data: 


1) 

w3.f 

(prepare  data  set  for  NADL) 

2) 

r2.f 

(prepeire  numerical  grid  coord  file) 

3) 

bilw.f 

(interpolate  onto  numerical  grid) 

4) 

wll.f 

(plot  selected  Wind  data) 

5) 

wil.f 

(plot  interpolated  Wind  field) 
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3.2.  The  3-D  Scalar  Field  Interpolation  Procedure 


The  Levitus  annual  mean  of  temperature  and  salinity  are  interpolated  on  the  following 
steps: 

1)  perform  a  2-D  extrapolation  of  the  Levitus  data  over  the  whole  NA  region  of  the  NADL. 
Extrapolation  is  done  using  the  ZGRID  algorithm  developed  by  J.  Taylor  et  al.  (1971). 

2)  perform  a  2-D  horizontal  bilinear  interpolation  of  T  emd  S  data  fields  at  each  z-level  of 
the  original  Levitus  data  set. 

3)  perform  a  1-D  vertical  interpolation  at  each  point  on  the  numericrd  grid.  The 

interpolation  is  done  using  a  cubic  spline  algorithm  from  Numerical  Recipes  (W.  H. 
Press  et  al.,1986). 

4)  compute  the  density  field  associated  with  the  T,  emd  S  data  distributions.  Density  is 
computed  from  the  UNESCO  ’80  equation  of  state  (Millero  and  Poisson,  1981). 

The  2-D  horizontal  interpolation  is  in  two  steps  because  of  the  coaurseness  of  the  Levitus 
data,  emd  the  misfit  between  the  numerical  grid  boundaries  and  the  geographiceJ  coastlines. 
The  extrapolation  over  the  NA  box  guarantees  that  not  many  undefined  data  points  are 
close  to  the  boundaries  of  the  numerical  grid  and  facilitates  application  of  the  bilinear 
interpolation  algorithm. 

The  separation  between  horizontal  and  vertical  interpolation  has  been  suggested  to  accel¬ 
erate  studies  of  the  sensitivity  of  the  model  to  the  vertical  grid  resolution  (i.e.,  we  expect 
to  conduct  several  experiments  configured  in  the  same  horizontal  grid  but  with  different 
vertical  structure).  The  2-D  interpolated  temperature,  salinity,  and  density  at  selected 
z-levels  are  depicted  in  Fig.  4. 

The  vertical  interpolation  algorithm  requires  the  definition  of  the  depth  at  each  numerical 
grid  point  and  can  be  applied  easily  to  any  general  vertical  coordinate  transformation 
adopted  by  the  ocean  circulation  model. 
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It  is  known  that  in  a  generalized  vertical  coordinate  system,  s,  conforming  to  the  bottom 
topography,  computation  of  the  pressure  gradient  forces  are  complicated  amd  might  intro¬ 
duce  sensible  errors  in  the  evolving  solution  (Haney,  1991).  In  order  to  improve  the  model 
performance,  it  is  customziry  to  decompose  the  density,  pxott  and  the  constituent  Tjot,  and 
Sxot  as  follows: 

PTot  =pbar(z)  +  p(^,rj,  s,t) 

Trot  =I^ar(z)  4-  T(^,7?,s,t) 

Sxot  =-^ar(z)  +  S($,?7,s,t) 

where  z  =  z(^,  t},s),  pbar,  Tbar,  and  Sbcir  represent  the  state  of  the  quiescient  ocean,  p,  T, 
and  S  the  perturbated  fields. 

To  augment  the  accuracy  of  the  density  and  pressure  forces  computations  and  increase 
the  stability  of  the  model,  it  is  necessary  to  ensure  that  the  perturbated  density  is  weakly 
stratified  with  respect  to  the  basic  state.  Since  the  problem  is  associated  with  the  numerical 
rather  than  the  physical  formulation  of  the  ocean  model,  the  basic  state  is  defined  ais  the 
mean  value  of  the  total  field  distributions  at  each  given  z-level.  (Since  the  equation  of  state 
is  not  linear,  pbar  csmnot  be  computed  from  Tbar  sind  Sbar,  directly).  Thus,  in  performing 
the  2-D  interpolation,  the  algorithm  also  computes  the  mean  values  for  the  defimtions  of 
the  basic  state.  Fig.  5.  illustrates  the  basic  state  of  the  quiescent  ocean  relative  to  the 
numerical  grid. 

The  numerical  algorithm  for  the  computation  of  the  density  field  also  includes  a  routine  to 
verify  that  the  density  field  is  stable  stratified,  and  eventually  performs  a  vertical  mixing 
of  the  constituent  fields. 

Five  programs  to  process  the  Levitus  data: 

1)  glll.f  (prepares  input  data  coordfile) 

2)  zgi.f  (extrapolation  of  input  data) 

3)  bilt.f  (interpolates  onto  numerical  grid) 
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4)  rhl.f  (computes  density  field) 

5)  vitl.f  (performs  vertica.’  interpolation) 

Five  programs  to  plot  Levitus  data: 

1)  tll.f  (plots  Levitus  data  on  the  NA  domeiin  at  a  given  z-level) 

2)  til.f  (plots  interpolated  2-D  data  at  a  given  z-level) 

3)  cil.f  (color  plots  of  2-D  interpolated  at  a  given  z-level  ) 

4)  dts.f  (plot  vertical  profile  of  the  basic  state) 

5)  tsl.f  (plots  T-S  diagram) 

4.  PLOTTING  FACILITIES 

Plotting  facilities  eire  for  the  display  and  presentation  of  the  data  sets  and  have  been 
described  in  the  previous  sections.  Plots  are  made  using  NCAR  Graphics  software.  The 
plotting  utilities  include  options  for  color  graphics  (not  presented  in  this  document)  and 
display  of  T-S  diagrams.  The  latter  facility  is  for  a  prelimineury  eveduation  of  the  water 
masses  present  in  the  data  sets,  to  trace  their  evolution  according  to  the  model  solution. 
Fig  6.  displays  the  T-S  diagram  of  the  upper  50m.  _ 

5.  CONCLUSIONS 

The  Hellerman’s  wind  and  Levitus  climatological  data  sets  were  acquired  and  processed 
for  preliminary  evaluations  of  an  eddy-resolving  oceein  circulation  model  for  the  NA. 

The  goal  is  to  implement  an  interpolation  package  that  can  be  applied  easily,  is  numerically 
efficient,  yet  preserves  the  dominant  featvires  of  the  original  records.  The  new  processed 
data  sets  are  part  of  the  initial  conditions  and  forcing  mechanisms  to  be  appled  in  numerical 
experiments  that  will  verify  the  physical  accuracy  and  numerical  efficiency  of  the  models. 
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FIGURES 

1)  The  North  Atlantic  dornmn  for  the  NADL,  and  the  orthogonal  boimd8iry  fitted  nu¬ 
merical  grid  used  in  this  document. 

2)  The  Hellerman  annual  wind  original  data:  (a)  The  longitude  (x)  component;  (b)  the 
latitude  (y)  component;  (c)  the  wind  stress  vector  directior.s  and  magnitude. 

3)  The  Hellerman  annual  wind  interpolated  data.  Same  as  Fig.  2. 

4)  The  2-D  interpolated  Levitus  data  at  the  z-levels:  50,  5C0,  1000  m.:  (a)  temperature; 
(b)  salinity,  and  (c)  density. 

5)  The  basic  state  of  the  quiescent  ocean  over  the  nvmerical  grid. 

6)  The  T-S  diagrams  of  the  upper  50m. 

7)  The  diagram  of  the  files  generated  by  the  interpolation  packages. 


8 


NORTH  ATLANTIC 
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Original  Wind  X  -  Annual  Mean 


CONTOUR  FflOH  -1.«000  tij  ’  ')000  CONTOUR  INTERVAL  OF  0.10000  PTI3  3I.  -0  5036 

Figure  2a 
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Origiiiiii  Wind  Y  Annual  Mean 
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Figure  2b 
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Wind  X  -  Annual  Mean 


CONTOUR  FROM  -.000005  TO  .000145  BY  .00001 
Figure  3a 
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CONTOUR  FROM  -.0000076  TO  .0000476  BY  .000006 
Figure  3b 
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Wind  with  magnitude-  Annual  Mean 
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Figure  3c 
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CONTOUR  FROM  0  TO  27  BY  1 

Figure  4a(l) 
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Temperature  -  level  13  500 


CONTOUR  FROM  2  TO  17  BY  1 

Figure  4a(2) 
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CONTOUR  FROM  3.6  TO  10.8  BY  .4 

Figure  4a(3) 
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Salinity  —  level  4  50 


CONTOUR  FROM  31.75  TO  37.25  BY  .25 

Figure  4b(l) 
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Salinity  -  level  13  500 


CONTOUR  FROM  34.6  TO  36.3  BY  .1 

Figure  4b (2) 
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CONTOUR  FROM  34.6  TO  36.2  BY  .1 

Figure  4b (3) 
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Density  -  level  4  50 


CONTOUR  FROM  23.5  TO  27,5  BY  .25 

Figure  4c (1) 


Density  -  level  13  500 

•/ 


CONTOUR  FROM  26.5  TO  26  BY  .05 

Figure  Ac (2) 
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CONTOUR  FROM  27.44  TO  27.76  BY  .02 


Figure  4c(3) 
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Wind  File  Diagram 


global  wind  data  set 


wind  on  the  NA  domain 


wind  interpolated  on  numerical  grid 


Levitus  File  Diagram 


Levitus  data  on  the  NA  domain 


Extrapolated  Levitus  data 


Levitus  data  interpolated  on 
numerical  grid  and  z-levels 


Data  interpolated  on  the  3D  numerical  grid 
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Figure  7 


APPENDIX;  USER’S  MANUAL 


This  section  describes  the  interpolation  numerical  procedures  as  well  as  what  the  user  has 
to  do  to  configure  the  prograims.  Changes  that  aire  the  user’s  responsibility  are  concentrated 
in  the  pcirameter  files  of  each  program.  Routines  that  are  included  in  the  procedures  are 
zdso  p2irt  of  the  NAUL.  Files  referred  to  as  fort.*  zu-e  created  from  the  grid  generation 
softweire.  The  diagraims  of  the  files  generated  by  the  progrcims  au'e  depicted  in  Fig.  7. 


Al:  THE  WIND  PROCESSING  PROGRAMS 


1.1.  PROGRAM:  w3.f 


The  program  writes  a  subset  of  the  Hellerman  winds  for  the  N'A.  The  selected  lon/lat  coor¬ 
dinates  are  written  to  a  file.  The  dimension  of  the  NA  domain  are  written  into  indata2.par 


Parameter  file:  w3.par 

NN 

LI  L2  Ml  M2 


wind.d 
wind.xy 
wind.  11 


number  of  data  sets  in  input  data  file 
selected  lon/lat  range  for  the  NA,  format  4(i4) 
Ll,L2:  Ion  range  (max  -179  to  179) 

Ml, M2:  lat  range  (max  -89  to  89) 
input  data  file 
NA  wind  data  file 
lon/lat  coord  of  NA  domain 


Program  expects  L\  <  L2,  Ml  <  M2.  Jf  1.1  >  L2  program  selects  Ll  to  179,  then  1  to 
L2.  If  Ml  >  M2,  program  changes  the  order  M1&M2,  Ml  <  M2.  The  parameter  file 
allows  the  parameters  to  be  changed  without  recompiling. 

Input  file:  wind.d  format  6el3.6 

This  file  contains  the  annual  mean  and  12  month  wind  data.  For  only  one  or  more  months, 
form  the  input  file  by  concatenating  from  hrjann-sm.d,  ml.d,  ...,  ml2.d,  and  change  the 
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number  of  data  sets  in  the  parameter  file  to  the  appropriate  number.  The  annual  meain 
counts  as  one  data  set,  as  do  each  of  the  months. 

Output  flies:  wind.xy  format  el3.6 

wind. 11  format  el3.6 

indata2.par  paxauneter  file 

1.2.  PROGRAM:  r2.f 

This  program  extracts  the  numerical  grid  lon/lat  coordinates  from  the  input  file. 
Parameter  flle:  grid2.par 

L  size  in  x  direction  of  numerical  grid 

M  size  in  y  direction  of  numericcd  grid 

Input  flle:  fort.62 

Output  flle:  spem.U  format  el3.6 

1.3.  PROGRAM:  bilw.f 

This  program  performs  bilinear  interpolation  of  the  wind  components  onto  the  numerical 
grid.  The  interpolated  vectors  then  need  to  be  adjusted  to  the  numerical  (^,  r})  orthogonal 
coordinate  axes. 

Parameter  flle:  grid.par  (cfr.  r2.f) 

Parameter  flle:  indata2.par 

NX  size  in  x  direction  of  NA  domain 

NY  size  in  y  direction  of  NA  domain 
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Parameter  file: 

bilw.par 

wind.xy 
wind. 11 
spem.ll 
wind.xt 

NA  wind  component  data  file 
NA  domain  lon/lat  data  file 
numerical  grid  lon/lat  data  file 
interpolated  output  data  file 

Input  flies: 

wind.xy 

wind.U 

spem.ll 
fort.  60 

(needed  for  axes  adjustment) 

Output  file: 

wind.xt 

format  12(el2.6) 

Routine: 

trgd2d 

(2-D  bilinear  interpolation) 

1.4.  PROGRAM;  wll.f 

This  program  creates  the  gmeta  plotfile  (contours  and/or  vectors)  for  the  NA  wind  data 
sets. 

Compile  using:  f77  -c  -03  -cg89  -o  wll.o  wll.f 

ncargf77  -cg89  -u  -o  wIplot.exe  wll.o 

Parameter  file:  indata2.par  (cfr.  bilw.f ) 


Parameter  file:  wll.par 


NN 

A  A  A  A  A 


wind.xy 

NNN 


number  of  months  selected 
plot  contour  options 
X  plots  X  contours 

Y  plots  Y  contours 

V  plots  vectors 

P  plots  vector  magnitude  &  direction 
input  data  file 

selected  months  (in  ascending  order) 
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Running  wlplot.exe  again  with  a  modified  par2uneter  file  wll.par  will  create  a  new  gmeta 
plotfile. 

Input  file;  wind.xy 

Output  file:  gmeta 

1.5.  PROGRAM;  wil.f 

This  program  creates  the  gmeta  plotfile  (contours  and/or  vectors)  for  the  wind  data  sets 
interpolated  on  the  numerical  grid. 

Compile  using:  f77  -c  -03  -cg89  -o  wil.o  wil.f 

ncargf77  -cg89  -u  -o  wiplot.exe  wil.o  plots. o  saxpy.o  grafs.o 

Parameter  flies:  grid2.par  (cfr.  r2.f) 

wll.par  (cfr.  wil.f) 

Running  wiplot.exe  again  with  a  modified  parameter  file  wil.par  will  create  a  new  gmeta 
plotfile. 

Input  flies:  wind.xt 

fort.  3 

Output  file:  gmeta 

Routines:  NCAR  graphics 

cpshift  (for  contours) 
velvet  (for  vectors) 
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A2:  THE  LEVITUS  DATA  PROCESSING  PROGRAMS 


2.1.  PROGRAM;  gll.f 


This  program  prepares  the  lon/lat  coord  file  for  the  Levitus  data  restricted  on  the  NA 
domain.  The  veJues  of  NX,  NY  (size  of  the  NA  domjiin)  are  written  to  indataS.pau-. 
Parameter  file:  gll-par 

Ll  L2  Ml  M2  selected  lon/lat  range,  format  4(i4) 

Ll,  L2;  Ion  range  (majc  -179  to  179) 

Ml,  M2:  lat  rauige  (max  -89  to  89) 
lev  1.11  lon/lat  coordfile  of  selected  data 


Output  files:  levl.ll 

indata3.pax 


format  el 3.6 
parameter  file 


2.2.  PROGRAM;  zgi.f 

This  program  performs  zgrid  interpolation  (on  reguleir  grid)  of  the  3-D  Levitus  tempera¬ 
ture/salinity  data. 

Parameter  file;  indata3.par 

NX  size  in  x  direction  of  Levitus  NA  domain 

NY  size  in  y  direction  of  Levitus  NA  domain 


Parameter  file;  zgi.par 


vtemp.dat 

vtemp.zdat 


input  data  file 
extrapolated  data  file 


Input  file: 
Output  file: 
Routine: 


vtemp.dat  or  vsal.dat  format  8fl0.2 

vtemp.zdat  or  vsal.zdat  format  flO.2 

zgrid  (2-D  interpolation) 
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2.3.  PROGRAM; 


bilt.f 


This  program  performs  horizontal  (bilinear)  2-D  interpolation  of  the  data  on  the  numerical 
grid. 

Parameter  flies:  indata3.par  (cfr.  tm.f) 


Parameter  file:  gridS.pax 


L 

M 


size  in  x  direction  of  numerical  grid 
size  in  y  direction  of  numericeil  grid 


Parameter  file:  bilt.pax 


vtemp.zdat 

levl.U 

spem.ll 

vtemp.zout 


input  data  file 

input  lon/lat  coordfile 

numerical  grid  lon/lat  coordfile  (cfr.  r2.f) 

interpolated  output  data  file 


Input  files:  vtemp.zdat  or  vsal.zdat  format  flO.2 

levl.U 
spem.ll 
fort.  3 

Output  flle:  vtemp.zout  or  vsal.zout  format  i3  el3.6  12(el2.6) 

level  meem  data 

Routine:  trgd2d  (2-D  bilinear  interpolation) 


2.4.  PROGRAM:  rhl.f 


This  program  computes  the  density  from  temperature  and  seJinity  2-D  interpolated  data. 
Parameter  flle:  grids. par  (cfr.  bilt.f) 


Parameter  flle: 

vtemp.zout 

vsal.zout 


rhl.par 

2-D  interpolated  temperature  data  file 
2-D  interpolated  saJinity  data  file 
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vdens.zout 


density  output  data  file 


Input  flies: 


vtemp.zout 

vsal.zout 


Output  file:  vdens.zout  format  i3  el3.6  12(el2.6) 

level  mean  data 


2.5.  PROGRAM:  vitl.f 

This  program  performs  vertical  (spline)  interpolation  of  the  3-D  data  on  the  numerical 
grid. 


Parameter  file:  grid3.par 


(cfr.  bilt.f) 


Parameter  flle:  vinterp.par 


XX 

XX 

N 


parameter  he 

pajzimeter  9 

number  of  levels  in  model 


Parameter  flle:  vitl.par 


vtemp.zout 

vitemp.zout 

vitbar.zout 


interpolated  (horizontal)  input  data  file 
vertically  interpolated  output  data  file 
vertically  interpolated  beisic  state  data  file 


Input  files: 


vtemp.zout,  vsal.zout,  or  vdens.zout 
fort. 3 


Output  flle:  vitemp.zout,  visal.zout,  or  videns.zout  format  12(el2.6) 

vitbar.zout,  visal.zout,or  vidbar.zout 


Routine: 


splint  (spline  interpolation) 


2.6.  PROGRAM:  tll.f 


This  program  creates  gnieta  plotfile  (contours)  for  the  Levitus  temperature  and  scilinity 
data  on  the  NA  domain. 

Compile  using:  f77  -c  -03  -cg89  -o  tll.o  tll.f 

ncairgf77  -cg89  -u  -o  tlplot.exe  tll.o 

Parameter  file:  indata3.paj  (cfr.  gll.f) 

Parameter  file:  til. pax 

NN  number  of  levels  selected 

A  T  for  temperature,  S  for  salinity 

vtemp.zdat  input  data  file 

NNNN  selected  levels  (ascending  order,  in  meters) 

Running  tlplot.exe  again  with  a  modified  parameter  file  til. par  will  create  a  new  gmeta 
plotfile. 

Input  file:  vtemp.zdat  or  vsal.zdat 

Output  file:  gmeta 

Routines:  NCAR  graphics 

cpcnrc 

2.7.  PROGRAM:  til.f 

This  program  creates  gmeta  plotfile  (contours)  for  data  sets  interpolated  on  the  numerical 
grid. 

Compile  using:  f77  -c  -03  -cg89  -o  til.o  til.f 

ncargf77  -cg89  -u  -o  tiplot.exe  til.o  plots.o  saxpy.o 
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Parameter  file: 


grids. par 


(cfr.  bilt.f) 


Parameter  file:  til. pax 


NN 

A 

vtemp.zout 

NNNN 


number  of  levels  selected 
T  for  temperature,  S  for  sedinity 
interpolated  input  data  file 
selected  levels  (ascending  order,  in  meters) 


Running  tiplot.exe  again  with  a  modified  parameter  file  til.peur  will  create  a  new  gmeta 
plotfile. 

Input  flies:  vtemp.zout  or  vszJ.zout 

fort. 3 


Output  file:  gmeta 

Routines:  NCAR  graphics 

cprect 


2.8.  PROGRAM:  cil.f 

Interactive  program,  creates  color  gmeta  plotfile  (contours)  for  the  data  interpolated  on 
the  numericzJ  grid. 

Compile  using:  f77  -c  -03  -cg89  -o  cil.o  cil.f 

ncargfr?  -cg89  -u  -o  ciplot.exe  cil.o  plots.o  saxpy.o 

Parameter  file:  grid3.par  (cfr.  bilt.f) 

Parameter  file:  cil.pjir 

NN  number  of  levels  selected 

NNNN  selected  levels  (ascending  order,  in  meters) 
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Running  ciplot.exe  ageiin  with  a  modified  parameter  file  cil.par  will  create  a  new  gmeta 
plotfile. 


Input  file: 

vtemp.zout,  vsal.zout,  or  vdens.zout 

Progrzim  prompts  for  following: 

A 

t  (temperature),  s  (scJinity),  or  r(density) 

A(30) 

data  filename 

A(26) 

title  of  plot 

NNN 

minimum  contour  level 

NNN 

maximum  contour  level 

Different  color  scales  are  used  for  temperature  aind  sa’inity: 
t:  grey/blue/green/red;  s,  or  r:  yellow /green/red/brown 

The  maximum/minimum  contour  levels  specify  the  reinge  of  contour  values  to  be  plotted. 
These  need  not  correspond  to  the  max/min  of  the  data  If  0  is  entered  for  both  max  and 
min  contour  levels,  then  the  program  uses  the  data  mzix  and  min.  Max/min  data  and  label 
values  for  the  contours  are  dumped  to  file  cil.out. 

Output  files:  gmeta 

cil.out 


2.9.  PROGRAM:  dts.f 

This  program  creates  a  ’gmeta’  plotfile  of  the  vertical  profile  of  temperature,  sedinity,  and 
density  basic  state  distribution. 

Compiled  using:  f77  -c  -03  -cg89  -o  dts.o  dts.f 

ncargf77  -cg89  -u  -o  dtsplot.exe  dts.o 

Parameter  file:  gridS.par  (cfr.  bilt.f) 
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Parameter  file: 


dts.paj 


temperature  data  file 
salinity  data  file 
density  data  file 


vtemp.zout 

vsal.zout 

vdens.zout 

Input  flies;  vtemp.zout 

vsal.zout 

vdens.zout 

Output  file;  gmeta 


2.10.  PROGRAM;  tsl.f 


This  program  creates  a  T-S  diagram  for  the  2-D  interpolated  data  sets  on  z-levels. 
Compile  using:  f77  -c  -03  -cg89  -o  tsl.o  tsl.f 

nceirgf??  -cg89  -u  -o  tsplot.exe  tsl.o 


Parameter  file;  grids. par 


(cfr.  bilt.f) 


Parameter  file;  gplot.par 


XL 

RX 

BY 

TY 


minimum  salinity  plotted 
maximum  salinity  plotted 
minimum  temperature  plotted 
maximum  temperature  plotted 


Parameter  flle;  tsl.pax 

NN 

NNNN 

NNNN 

vtemp.ozut 

VSEll.ZOUt 


interval  between  points  plotted 
if  1,  every  point  in  data  files  is  plotted 
minimum  depth  (in  meters)  of  point  plotted 
maximum  depth  (in  meters)  of  points  plotted 
interpolated  temperature  data 
interpolated  salinity  data 


Input  files; 


vtemp.zout 

vsal.zout 


Output  file;  gmeta 
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